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ABSTRACT 

In galaxy clusters the equilibria of the intracluster plasma (ICP) and of the gravitationally dominant dark 
matter (DM) are governed by the hydrostatic and the Jeans equation, respectively; in either case gravity is 
withstood by the corresponding, entropy-modulated pressure. Jeans, with the DM 'entropy' set to K oc r a 
and a « 1.25-1.3 applying from groups to rich clusters, yields our radial a-profiles; these, compared to 
the empirical NFW distribution, are flatter at the center and steeper in the outskirts as required by recent 
gravitational lensing data. In the ICP, on the other hand, the entropy run k(r) is mainly shaped by shocks, as 
steadily set by supersonic accretion of gas at the cluster boundary, and intermittently driven from the center 
by merging events or by active galactic nuclei (AGNs); the resulting equilibrium is described by the exact 
yet simple formalism constituting our ICP Supermodel. With a few parameters, this accurately represents 
the runs of density n{r) and temperature T(r) as required by up-to-date X-ray data on surface brightness and 
spectroscopy for both cool core (CC) and non cool core (NCC) clusters; the former are marked by a middle 
temperature peak, whose location is predicted from rich clusters to groups. The Supermodel inversely links 
the inner runs of n(r) and T(r), and highlights their central scaling with entropy n c oc k~ l and T c oc £|?- 35 , to 
yield radiative cooling times t c w 0.3(fc c /15keV cm 2 ) 12 Gyr. We discuss the stability of the central values so 
focused: against radiative erosion of k c in the cool dense conditions of CC clusters, that triggers recurrent AGN 
activities resetting it back; or against energy inputs from AGNs and mergers whose effects are saturated by the 
hot central conditions of NCC clusters. From the Supermodel we derive as limiting cases the classic polytropic 
/3-models, and the 'mirror' model with T(r) oc <r 2 (r) suitable for NCC and CC clusters, respectively; these 
limiting cases highlight how the ICP temperature T(r) strives to mirror the DM velocity dispersion er 2 (r) away 
from energy and entropy injections. Finally, we discuss how the Supermodel connects information derived 
from X-ray and gravitational lensing observations. 

Subject headings: Dark matter — galaxies: clusters: general — gravitational lensing — methods: analytical 
— X-rays: galaxies: clusters 



1. INTRODUCTION 

Diffuse baryons and dark matter (DM) constitute the major 
components of clusters and groups of galaxies, with the for- 
mer energized and shining by continual struggle against the 
latter. 

The DM accounts for some 6/7 of the total masses M ~ 
10 13 - 10 15 M Q from poor groups to rich clusters, making 
for average densities p ~ 10~ 26 g cm" 3 with its constituent 
('collisionless') particles entertaining little or no interactions 
other than gravity. Thus the DM sets the overall gravitational 
wells virialized within radii R up to a few Mpcs, where all 
bodies in dynamical equilibrium - from single particles to 
whole galaxies - possess or acquire an 1-D velocity disper- 
sion a 2 ks GM/5R ~ 10 3 km s _1 . 

The bulk of the baryons, to a fraction again close to 6/7, 
is found in the diffuse form of a hot intracluster plasma (the 
ICP), mostly comprised of protons with the neutralizing elec- 
trons at number densities n ~ 10~ 3 cm -3 , and equilibrium 
temperatures IcbT rs m p a 2 /2 ~ 5 keV (Icb being the Boltz- 
mann constant and m p the proton mass) well above most 
ionization potentials. This we know since Cavaliere et al. 
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(1971) found in the first Uhuru data clear evidence of a new 
class of bright X-ray extragalactic sources associated with the 
deep, stable gravitational wells of the galaxy systems, emit- 
ting from the ICP they contain thermal bremsstrahlung pow- 
ers L x ~ 2 x 1 0" 27 n 2 R 3 T 1 / 2 ~ 1 42 — 1 45 erg s" 1 . The notion 
has been nailed down beyond all doubts by observations of the 
extended nature of these sources (Gursky et al. 1972) and of 
high excitation, coronal-like lines (Mitchell et al. 1976; Ser- 
lemitsos et al. 1977; see also Sarazin 1988) that also pointed 
toward definite, somewhat subsolar metallicities. 

Such temperatures and densities make the ICP an extremely 
good plasma, in fact the best in the Universe ever, as its con- 
stituent particles in the DM gravitational wells acquire a large 
kinetic relative to electrostatic energy (at mean separations d = 
n" 1 / 3 ~ 10 cm), their ratio being of order ksT /e 2 n^ 3 ~ 10 12 ; 
this astounding value is to be compared with its counterparts: 
10 3 in stellar interiors, or 3 x 10 s in the pre-recombination 
Universe. This holds despite gravity being so exceedingly 
feeble at microscopic levels as to attain a mere Gm 2 /e 2 ~ 

8 x 10~ 37 of the strength that marks the electromagnetic inter- 
actions; it holds because the plasma condition 

10 12 - k B T/e 2 n l/3 = Gm 2 p /e 2 x d/lOR x N (1) 

is dominated by the huge number M = M/m p ~ 10 73 express- 
ing in proton units the total DM mass with its overwhelm- 
ing gravity. As a result, the ICP at microscopic scales con- 
stitutes a very simple, nearly perfect gas of particles with 3 
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degrees of freedom and effective mass pm p with p, « 0.6. At 
intermediate scales of some 10 kpc these sense mostly the 
electromagnetic interactions; by the latter the electrons ab- 
sorb or emit radiation, while the protons share their energy 
over a mean free path \ pp w 10(fcBr/5keV) 2 (n/T0~ 3 cm~ 3 ) _1 
kpc, with electrons following suit over some 40X pp toward 
full thermal equilibrium. 

By the same token, the macrophysics at cluster scales gets 
interestingly complex, as the ICP constitutes a faithful archive 
preserving memory of the different energy inputs occurring 
throughout sizes of Mpcs and over timescales of 10 Gyr. We 
shall see that the ICP physics is governed, in a nutshell, by the 
interplay of electromagnetic interactions transferring energy 
to the plasma, and of the DM bulk gravity proceeding to blend 
and readjust it over cluster scales. 

To disentangle these processes, it will prove technically 
convenient to use two synthetic and formally analogous quan- 
tities. As for the ICP, the adiabat k = ksT /n 2 / 3 is straight- 
forwardly related to the thermodynamic specific entropy s = 
3£fi/2 x log£+ const. Usually named 'entropy' for short, this 
quantity is endowed with time-honored, effective properties 
like: increasing (decreasing) by energy gains (losses) other 
than adiabatic compressions (expansions); and controlling the 
ICP settlement into gravitational wells. 

As to the DM, the analogous quantity K = a 2 /p 2 / 3 is in- 
creasingly found to play - in spite of all the traditional objec- 
tions leveled to defining 'entropy' in a collisionless medium 
dominated by long range self-gravity like the DM - similar 
roles to a true entropy; in fact, it increases during the fast col- 
lapse with associated major mergers that set up the halo bulk, 
and stays put during the subsequent slow mass accretion pro- 
ducing a quiet development of the outskirts from the inside 
out. These two stages have been recently recognized in so- 
phisticated A^-body simulations that follow the development 
of halos embedding galaxies or galaxy systems from initial 
cold DM perturbations (e.g., Zhao et al. 2003; Diemand et al. 
2007). From simulations, during quiet stages the simple pow- 
erlaw run K(r) oc r a is found to apply remarkably well with 
definite slopes close to a « 1.25 throughout the structure's 
main body, as first stressed by Taylor & Navarro (2001) and 
confirmed by many others (e.g., Dehnen & McLaughlin 2005; 
Hoffman et al. 2007; Ascasibar & Gottlober 2008; Vass et al. 
2009; Navarro et al. 2009). 

We shall find these two entropies to be most useful in com- 
puting and relating the macroscopic, static equilibria of DM 
and baryons in the same gravitational wells as set by the for- 
mer. This is perceived on just inspecting their related equilib- 
rium conditions for the density runs n(r) and p(r), namely, the 
hydrostatic vs. Jeans equation in the form 
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pm p n dr 



(n 5/3 /fc) = - 



GM(< r) _ 1 d 
r 2 p dr 



(P 5/3 K) . (2) 



While the rich contents of the second equality have been dis- 
cussed in detail by Lapi & Cavaliere (2009a [hereafter LC09], 
see also § 2 for a recap), in the present paper we will focus on 
the first equation. To this purpose, we make use of a physical 
model for k(r) to derive the detailed pressure structure of the 
ICP. 

Our approach develops the line extensively pursued in 
the literature at increasingly sophisticated levels. A simple 
isothermal or polytropic state has been adopted by, e.g., Cav- 
aliere & Fusco-Femiano (1976, 1978), Balogh et al. (1999), 
Dos Santos & Dore (2002), and Ostriker et al. (2005). Re- 



fined models were based on a more articulated entropy run, 
consisting of an outer ramp produced by gravitational accre- 
tion shocks going on in the aftermath of cluster formation (see 
Tozzi & Norman 2001, Voit et al. 2003, Lapi et al. 2005 
[hereafter LCM05]), plus a central entropy floor of nongravi- 
tational origin (e.g., Babul et al. 2002; Voit et al. 2002; Voit 
2005). The latter may be contributed by a number of pro- 
cesses: cooling that selectively removes some low-entropy 
gas from the inner regions (see Bryan 2000; Voit & Bryan 
2001; McCarthy et al. 2004); preheating of the ICP due to su- 
pernovae (SNe)/active galactic nuclei (AGNs) within the clus- 
ter progenitors (see Cavaliere et al. 1997; LCM05; McCarthy 
et al. 2008); the direct action within the cluster of central 
AGNs (Valageas & Silk 1999; Wu et al. 2000; Scannapieco 
& Oh 2004; LCM05). Additional gravitational events in the 
form of major mergers (Balogh et al. 2007; McCarthy et al. 
2007) that now and then punctuate the equilibrium state of 
formed clusters (Cavaliere et al. 1999) may reach down to the 
center and enhance the entropy there. 

Here we improve over these previous studies in the fol- 
lowing respects. As to the DM potential well, we adopt the 
physical a-profiles obtained from solving the Jeans equation 
rather than assuming the empirical NFW fit (see § 2). As to 
the ICP equilibrium, we present a novel analytical solution 
for the ICP density and temperature runs in terms of the ICP 
entropy distribution (see § 3). The distribution that we use 
comprises an outer ramp and a central floor, marked by two 
parameters physically assessed in terms of steady, self-similar 
accretion shocks plus additional energy inputs due to central 
AGNs and mergers (see § 3.1). We insert the entropy distri- 
bution in our equilibrium solution to obtain the ICP 'Super- 
model', that effectively represents with two ICP parameters 
the extended profiles of the X-ray surface brightness and tem- 
perature for either classes: the cool core (CC) and the non cool 
core (NCC) clusters, see § 4. We discuss the stability against 
cooling and feedback of the ensuing central conditions (see 
§ 5). We recover classic models for the ICP distribution of 
density and temperature as limiting cases of the Supermodel 
valid in different radial ranges and for different central en- 
tropy levels (see § 6). We provide specific predictions for 
future X-ray observations of the ICP, and model-independent 
tests of the entropy runs underlying the ICP Supermodel (see 
§7). 

Throughout this work we adopt a standard, flat cosmol- 
ogy (see Spergel et al. 2007) with normalized matter density 
VLm = 0.27, dark energy density f^A = 0.73, and Hubble con- 
stant H = 72 km s" 1 Mpc" 1 . 

2. a-PROFILES FOR THE DM HALOS 

From LC09 we recall for use in the present paper a number 
of basic features in the equilibrium of the DM halos embed- 
ding the ICP, as described by the Jeans equation with entropy 
K oc r a (see § 1); in compact form this reads 
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(3) 



With a = dlog/f/dlogr set to a constant, this yields the 
changing density slope -f(r) = -dlogp/dlogr in terms of the 
increasing ratio to a 2 (r) of the halo circular velocity squared 
v 2 (r) = GM(< r)/r that provides normalization and running 
estimates for the gravitational potential. The above equation 
provides the detailed cluster gravitational well, within which 
the ICP is to adjust. 
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FIG. 1. — Top panel: radial runs of the DM density (solid) and entropy 
(dotted); bottom panel: radial runs of the DM velocity dispersion (solid) and 
circular velocity (dotted). Both panels refer to the a-profiles with a = 1.27 
(green lines) and a = 1.3 (red lines); the profiles are normalized to unity at 
the radius ro where the density slope 70 applies. The behaviors in the in- 
ner regions, main body and outskirts are highlighted by black dashed lines. 
The radii r_i where the density slopes equal -2 (used in defining the con- 
centration c, see § 2), r m where the velocity dispersion peak, and r p where 
the circular velocity peak are highlighted as squares, filled and empty circles, 
respectively. 

The physically relevant values of a are pinned down from 
the structures' cosmogonic development. Based on the 
straightforward scaling laws a 1 oc M/r and r oc M/M 2 / 3 , it 
is seen that the entropy is to scale as K oc rM 1 / 3 , a run clearly 
implying a radial slope a > 1. The detailed time depen- 
dence of the accretion rate M/M and the related quantities 
may be derived semianalytically within the standard ACDM 
cosmogony; the result for the allowed range is a = 1 .25 - 1 .3 
from groups to massive clusters, narrowed down to 

aw 1.27-1.3 (4) 

for average masses M ~ 10 14 -10 15 Mq from poor to rich clus- 
ters of main interest here. In the structure's development these 
values obtain at the transition epoch from the stage of fast col- 
lapse to that of slow accretion (relative to the running Hub- 
ble time), when the potential well attains its maximal depth 
marked by Mt/M ~ 1. Such a two-stage development turns 



out to be in tune with the intensive, detailed A^-body simula- 
tions recalled in § 1. In addition, the simulations support a 
value of a in the range above, which stays put from the center 
throughout the halo bulk as the late quasi-equilibrium config- 
uration develops from the inside out. 

These values may be inserted into the Jeans equation, to 
find (as pioneered by Taylor & Navarro 2001 and Dehnen 
& McLaughlin 2005) all the viable solutions for p(r), that 
we name 'a-profiles' and illustrate in Fig. 1. These fea- 
ture a monotonic radial run satisfying physical central and 
outer boundary conditions, i.e., zero gravitational force (cor- 
responding to a flat minimum of the potential) and finite 
(hence definite) overall mass, respectively. 

The ensuing behavior of the a-profiles, basic to the ICP 
equilibrium, is illustrated in the top panel of Fig. 1, and high- 
lighted by the analytic expressions of the slopes 

3 3 
7«=5«. 7o = 6-3a, 7& = -(l + a). (5) 

These start from the central (r — > 0) value j a « 0.76-0.78, 
progressively steepen to 70 « 2.19-2.1 at the point ro that 
marks the halo main body, and steepen further into the out- 
skirts to the value 72, « 3.41-3.44 at around the virial ra- 
dius R before going into a final cutoff. The a-profiles are 
seen (see LC09) to correspond to a maximal value n a -it(ct) = 
v 2 /tx 2 w 2.6-2.5 for the relative gravitational pull at the 
point r = r p > ro where v 2 (r) peaks (see also bottom panel of 
Fig. 1). After Eq. (3) this also implies at r p a maximal slope 
7 P = 3(a+fc cri t)/5« 2.32-2.28. 

At variance with the empirical NFW formula (Navarro et 
al. 1997) that features angled central potential or pressure and 
diverging mass, for the a-profiles the inner slope is consider- 
ably flatter, and the outer one is steeper as to result in a defi- 
nite overall mass. The radial range r > r_2 where the density 
profile is steeper than a reference slope 7 = 2 may be specified 
in terms of the usual concentration parameter c = R/r_2, that 
may be viewed as a measure of central condensation and/or 
outskirts' extension. In up-to-date numerical simulations (see 
Zhao et al. 2003; Diemand et al. 2007) this is found to take 
on value c w 3.5 at the transition redshift z,, and to increase 
thereafter to current values c w 3.5 (1+z,) up to c « 10 for 
the fraction about 10% of rich clusters with early transition 
epoch Zt ~ 1-5. Instead, values c k 4-5 apply to the much 
more numerous clusters with recent transition epochs. 

Another basic feature of the a-profiles is provided by the 
peaked run of the gravitationally acquired dispersion 

a\r)KK{r) P 2 l\r), (6) 

see bottom panel of Fig. 1. Towering above the central and 
the outer drops related to the cold nature of the DM, the peak 
results from the central rise of K(r) and the outer steep falloff 
of p(r), and will have a striking counterpart in the ICR The 
characteristic values of a 2 for the a-profiles are set in terms 
of minimal values for a 1 = v 2 /K C iit', this is related to limited 
randomization during accretion of the infall kinetic energy 
gained by the initially cold DM particles. 

These a-profiles actually depend weakly on a in the named 
range (see Fig. 1), and though derived from the isotropic 
Jeans Eq. (3) prove to be stable against addition of reasonable 
anisotropies. The latter are described by the standard Bin- 
ney (1978) parameter (3, which numerical simulations suggest 
to increase outwards from central values /3(0) > -0.1 mean- 
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ing weakly tangential anisotropy, toward /3 < 0.5 that implies 
prevailing radial motions (see Dehnen & McLaughlin 2005; 
Hansen & Moore 2006; also H0st et al. 2009). Such den- 
sity profiles turn out to be slightly flattened at the center and 
considerably steepened into the outskirts. 

The a-profiles are also stable against educated variations 
in the condition a = const. Actually, the cosmogonic buildup 
given by our semianalytic computation yields a slow decrease 
of a(r) into the outskirts, as these develop from the inside out 
at late times after the transition. Again in keeping with the 
simulations, the outer development does not affect the inner 
gravitational potential and the equilibrium described by the 
Jeans equation, that (with its central boundary condition) also 
works from the inside out. In fact, for the present use we 
have computed the density profiles associated with a(r), and 
checked these to be very close to the parent a-profiles but 
for being somewhat steeper into the outskirts, with densities 
lower by about 15% in the vicinity of the virial radius. 

It is to be stressed that the a-profiles for M ~ 10 14 — 
10 15 M Q , especially with full a(r) and anisotropics, turn out 
to provide an optimal fit to the surface density runs as de- 
rived from gravitational lensing observations in around mas- 
sive clusters that just require central slopes flatter and outer 
slopes steeper than the standard NFW (cf. Broadhurst at al. 
2008). Such a fit for A1689 requires high values c ~ 10 and 
early transition epochs z t w 1.5, related conditions that the 
two-stage development predicts to occur in about 10% of the 
rich clusters (Lapi & Cavaliere 2009b). 

We will adopt these a-profiles depending on the two DM 
key parameters a and c to describe the gravitational potential 
containing the ICP, to which we now turn. 

3. THE ICP EQUILIBRIUM 

We will focus on clusters in equilibrium conditions, in be- 
tween the punctuating violent mergers for which we refer 
the reader to the review by Markevitch & Vikhlinin (2007). 
When the DM halo is close to equilibrium, even more so 
will be the pervading ICP; in fact, the sound crossing time 
R/(5 ksT /3 [imp) 1 ! 1 is seen to be somewhat shorter than the 
dynamical time R/a on recalling from § 1 that kgT w pm p a 1 
holds. In these conditions the ICP is governed by the hydro- 
static, Jeans-like equation provided by the first and the second 
member of Eq. (2). In a compact form similar to Eq. (3) this 
writes 

3 3 

g=-a+-b, (7) 

once the ICP entropy run k(r) is given and used to factor 
out the product d(kn 5 ^)/dr in the original Eq. (2). In fact, 
the entropy slope a = dlogfe/dlogr enters Eq. (7) to yield 
the density slope g(r) = -dlogn/dlogr, concurring with the 
gravitational pull measured by the (squared) ratio 3b(r)/5 = 
3 pm p v 2 c {r)/5kBT{r) of the DM circular velocity (discussed in 
§ 2) to the sound speed. 

The above constitutes just a first order differential equation, 
that transforms to linear in terms of the variable « 2 / 3 (r). By 
textbook recipes (cf. Dwight 1961) this is amenable to a sim- 
ple solution in form of a quadrature, that is conveniently writ- 
ten as 

2 f 1 dr' 

T(r) = k(r)n 2/3 (r) = k 3 ' 5 (r)[l + -b R I — v 2 (r')r 3/ V)] . 

(8) 
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FIG. 2. — Top panel: the integral entering Eq. (8) is plotted as a function of 
r for the values k c = 5 X 1(T 2 (red), 2.5 X 1(T 2 (cyan), and (blue) used in the 
next Figures; the dotted lines represent the analytic fit provided in Appendix 
A. Bottom panel: The value of the integral at the center is plotted as a function 

--1/4 

of k c (solid line), and is compared with the powerlaw approximation k c 
(dotted line); the latter holds to better than 5% in the most relevant range 
kc?z 0.01-0.15, see §4. 

Here barred variables are normalized to their boundary value 
at r = R where b(r) takes on the value b R , while v 2 c {r) is taken 
from the a-profiles with its weak dependence on a. 

In Fig. 2 (top panel) we plot the integral appearing in the 
above equation, along with the simple analytical fit presented 
in Appendix A. Note three circumstances. First, 2b R /5 ks 1 
holds (e.g., 2b R /5 « 1.06-1.02 on using b R = 2.65-2.55, 
see Eq. [10]). Second, at the center the integral dominates in 
the square bracket; thus the memory of the boundary condi- 
tion is swamped while the integral is numerically found to 
scale as kc with the central entropy k c , see Fig. 2 (bot- 
tom panel). Third, on approaching the center where the r- 
dependence of the bracket is already saturating, Eq. (8) im- 
plies T(r) oc «"'(r) oc £ 3 / 5 (r). The issue to stress is that, rather 
than a model, the above Eq. (8) has the standing of a theorem 
in hydrostatics valid for the ICP of clusters close to equilib- 
rium. 

3.1. ICP entropy 

What actually needs physical modeling is the entropy run 
k(r). We base upon the notions that entropy is erased by 
radiative cooling on the timescale t c « 65(ksT /5keW)^ 2 
(«/10~ 3 cirf 3 ) -1 Gyr (Sarazin 1988), while substantial raising 
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requires shocks (see Appendix B) as are driven by supersonic 
outflows from center and set by inflows across the boundary. 

3.1.1. Entropy deposited at the boundary 

At the outer end, the mathematics of Eqs. (7) requires one 
boundary condition fixing b R . On the observational side, T(r) 
is found to decline slowly toward the virial radius (Molendi & 
Pizzolato 2001), and at such a rate it would take tens of Mpcs 
to decline smoothly fromkeV values to those some 10~ 2 times 
lower as prevailing in the external medium. So at r « R a dis- 
continuity is to occur and terminate the hot ICR at variance 
with the smooth if steeper decline of the DM density. In fact, 
most of the transition takes place across a few mean free paths 
\ pp at the accretion shocks produced when external gas super- 
sonically falls into the DM potential. 

Numerical simulations (e.g., Tormen et al. 2004; see their 
Fig. 8) show that accretion of the smooth gas and minor lumps 
making up a major fraction of the accreted mass drives a com- 
plex patchwork of shocks mostly comprised within an outer 
layer with thickness of order 10 _1 R; across such a layer the 
standard conservation laws of mass, momentum and energy 
may be applied leading to the classic Rankine-Hugoniot jump 
conditions; to within 10% the results are similar to a coher- 
ent accretion shock, roughly spherical and located at r ~ R 
(see LCM05 for a complete treatment, and Ettori & Fabian 
1998 for possibly delayed electron equilibrium). The impact 
of major lumps reaching down to the center will be dealt with 
in §3.1.3. 

The jump conditions take on a particularly simple form for 
cold inflow into rich clusters that cause strong shocks with 
Mach numbers squared A4 2 considerably exceeding unity 
(see Appendix B). In strong shocks lingering at r « R (and ex- 
panding along with the cluster's development), maximal con- 
version of infall energy occurs over a radial range of order 
X pp , to yield 

2 2 

k B T R = - ixm p v R A0tR , (9) 

where 2v\ A</> t R is the kinetic energy per unit mass of the gas 
freely falling across the potential drop down to r = R from the 
turning point where infall starts (see LCM05). Thus we find 



b R = 



(10) 



with values 2.65-2.55 corresponding to a = 1.27- 1.3 and 
A4> tR w 0.57-0.59. 

At the boundary the equilibrium Eq. (7) yields g = 3 (a + 
b R )/5. Meanwhile, powerlaw approximations k oc r a , n oc r~ g 
apply in the vicinity of r = R, so m(r) oc r 3 ~ g describes the ICP 
mass m in the outer layer, and correspondingly 



k OC m a / (3 - ?) OC m 5 «/3(5-a-fo S ) 



(11) 



obtains. 

On the other hand, pursuing the scaling for K(r) recalled 
in § 2, LC09 show that the DM entropy behaves as K oc 
M 4 / 3 /M 2 / 3 oc M 3 / 2 on considering that M oc f 4 / 5 holds in the 
standard ACDM cosmogony for z < 0.5 during the slow ac- 
cretion stage. As external gas and DM are accreted in cosmic 
proportion, in the outer layer m oc M applies, and it follows 
that the ICP entropy at the boundary can be expressed as 



Then on equating the exponents in Eqs. (11) and (12) one 
finds 

45-9fr g 27 + 6b R 

or= — — — and g R = — — — . (13) 
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With the values of b R discussed above, these yield the slopes 
a R w 1.1 and g R m2.2, both considerably flatter than the cor- 
responding DM values. Note from Eq. (13a) how the slope a R 
at the boundary decreases if b R is increased. We expect such 
a decrease to take place for clusters with high concentrations 
corresponding to shallow outer potentials (see Lapi & Cava- 
liere 2009b), that decrease the outer drops A^tR entering b R 
through Eq. (10). E.g., a value c « 10 (holding for a cluster 
with an early transition, see § 2) in place of the usual c « 4 
(holding for clusters with a recent transition) implies A(j) tR to 
lower to 0.47, b R to grow to 3.2, and a R to decrease to 0.85 
as may be the case for A1689 (see Lemze et al. 2008; Lapi 
& Cavaliere 2009b); meanwhile, from Eq. (13b) g R increases 
to 2.4. A flat entropy slope may be also produced when the 
boundary shock is weakened by substantial preheating of the 
infalling gas; such a condition is relevant for poor clusters and 
groups with comparatively low potentials (see LCM05), e.g., 
a w 0.8 applies to systems with M < 10 14 M Q with external 
preheating at levels around 1 /2 keV per particle. 

On the other hand, after Eqs. (13) high values of a with 
an upper bound at 45/19 w 2.4 (in keeping with the variance 
observed by Cavagnolo et al. 2009) correspond to low val- 
ues of b R , hence to large A0 t R; these imply flat densities and 
temperatures sustained to high values in the outskirts. 

3.1.2. Entropy stratification 

No other major sources or sinks of energy and entropy oc- 
cur inward of the boundary at a few Mpcs down to the central 
10 2 kpc. So throughout the ICP bulk, the entropy slope is 
to stay at its boundary value a = a R , and in rich clusters with 
standard concentration 



k(r) oc r 



l.i 



(14) 



k oc m 3 / 2 . 



(12) 



is to hold, a result similar to Tozzi & Norman (2001), LCM05, 
and Voit (2005). 

On the other hand, while the simple thermodynamics of 
the ICP maintains the run of k(r) in this powerlaw form, 
its more complex adiabatic readjustments within the gravita- 
tional well cause the density slope g(r) to flatten out inward; 
this is granted by Eq. (7), as the inequality b(r) < b R (i.e., 
v 2 (r) / 'ksT '(f) < 1) progressively strengthens away from the 
boundary. As a result, the ICP slope g(r) not only starts out 
but also stays generally flatter than DM's j(r) inward of the 
boundary, as given by 

g(r) = 3 [a + b(r)]/5 < 3 [a + K (r)]/5 = 7 (r) . (15) 

This implies the ICP density run n(r) at the boundary to paral- 
lel the DM's p(r) at some other point well inside the cluster's 
main body. 

3.1.3. Entropy in the central region 

Next we discuss the central value k c of the entropy and its 
origins. Recall from § 1 that the central entropy may be pro- 
duced by shocks driven by substantial merging events reach- 
ing the center (McCarthy et al. 2007; Balogh et al. 2007), 
and by AGNs residing in the central bulge-dominated galax- 
ies (see Valageas & Silk 1999; Wu et al. 2000; Cavaliere et 
al. 2002; Scannapieco & Oh 2004; LCM05). In addition, a 
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FIG. 3. — Radial runs of the ICP density, mass, entropy, temperature, cool- 
ing time, and polytropic index given by the Supermodel. Solid lines are for 
different values of the central entropy k c = 5 X 10 r 2 (red), 2.5 X 1CT 2 (cyan), 
and (blue); dotted lines refer to the 'mirror' model discussed in § 6.1 with 
fi = 0.75; dashed lines illustrate the underlying DM distribution taken from 
the a-profile with a = 1.27 and an average concentration value c = = 5 

(see § 2). In the panel representing T(r), we plot the DM velocity dispersion 
cr 2 (r) as a thick dashed line, and the circular velocity v 2 (r) as a thin dashed 
line; we also show average data for the CC (circles), NCC (squares), and 
UNC (triangles) classes from Leccardi & Molendi (2008). The hatched area 
outlines the central range hardly accessible to current resolutions. 



basal entropy level may be left over by SNe and AGNs that 
preheated the gas in the volume due to collapse into, or to 
accrete onto the cluster (Cavaliere et al. 1997, LCM05, Mc- 
Carthy et al. 2008). 

To begin with, consider preheating; this may be expressed 
in terms of entropy advected across the shock by the currently 
infalling gas (see Appendix B). However, such a process is 
known to be effective only within limits: if always strong, it 
would flatten the bulk entropy slope of most clusters to values 
a < 1 (see § 3.1.1 and Fig. 3), and produce really fiat entropy 
profiles in most groups at variance with the observations (see 
Balogh et al. 1999; Pratt & Arnaud 2003; Rasmussen & Pon- 
man 2004); unless tailored to cluster precollapse size, it would 
involve larger volumes and require high total energetics; fi- 
nally, the presence of Lya absorbers limits actual preheating 
to optimal redshifts z ~ 2 - 3 when it may catch the peak of 
star/AGN formation while still operating at moderate cosmo- 
logical densities (see Voit 2005). 

Consider then the shocks produced at cluster centers by ma- 
jor mergers. A main progenitor with mass already close to the 
present value M undergoes only a few more, if any merging 
shocks; each of these, in view of the considerable mass import 
and high infall velocity, can deliver large energies up to 10 64 
ergs or several keV per particle, and may produce levels of k c 



FIG. 4. — ICP projected distributions of emission-weighted temperature, 
spectroscopic-like temperature, X-ray surface brightness, and l.o.s. SZ effect. 
Linestyles as in Fig. 3. 

up to several 10 2 keV cm 2 (see Appendix B). With timescales 
for free radiative cooling of order 

t c « 0.3 ( K . ) ( ^-Y Gyr , (16) 
Vl5keVcm 2 / V5keVy 

(Voit & Donahue 2005), the above levels take several Gyrs to 
decrease substantially. At the other end, a much lower entropy 
state would rapidly run down to values lower still. So a key 
point is to secure intermediate but persisting entropy levels. 

Values of k c in the range 10-30 keV cm 2 as frequently 
observed despite fast radiative erasure point toward entropy 
inputs frequently refreshed; a promising path is provided by 
sizeable and recurrent energy injections by central AGNs, 
triggered into activity by renewed fueling in a low-entropy, 
dense environment of their underlying supermassive black 
holes (BHs). These are observed to undergo outbursts able 
to inject, nearly independent of cluster mass, energies AE up 
to 10 62 erg or a few keV per particle, yielding k c of several 
tens keV cm 2 . Cooling on the scale of 10 _I Gyr may be offset 
as the several supermassive BHs inhabiting the many bulge- 
dominated galaxies in the central region of a rich cluster alter- 
natively kindle up, and collectively recur over such timescales 
(see Nusser et al. 2006; Ciotti & Ostriker 2007; Conway & 
Ostriker 2008). In particular, outbursts lasting for a few 10 8 
yr can continuously drive out to 300 kpc pressurized blast- 
waves terminating into a shock with Mach numbers sustained 
at M 2 > 3, as computed in detail by LCM05 and observed 
by Forman et al. (2005), Nulsen et al. (2005), McNamara & 
Nulsen (2007), see also Appendix B. 

All these processes contribute to raise the entropy over an 
extended region around the center, adding to the outer run 
given by Eq. (14); the combined entropy profile may be de- 
scribed by the simple parametric expression (see Voit 2005 
and references therein) 

k(F) = ~kc + d-~kc)r a (17) 

represented in Fig. 3. Having computed the outer powerlaw 
slope a around 1.1 (see § 3.1.1) and having discussed central 
entropy levels k c from a few tens to several hundreds keV cm 2 , 
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FIG. 5. — X-ray surface brightness and temperature of the NCC cluster 
A2218. Circles represent XMM-Newton data (Zhang et al. 2008), while 
squares refer to preliminary Suzaku observations. For the sake of clarity, in 
the top panel we do not report the formal error bars; these are generally small, 
except for the outermost data limited by sensitivity, and for the innermost 
ones limited by resolution. The solid line represents the Supermodel for a 
given DM a = 1.27-profile with concentration c = 3.5, and ICP entropy with 
slope a=l.l and central value k c = 10 _1 . 



we next use these two ICP key parameters in our equilibrium 
'Supermodel'. 

4. THE SUPERMODEL FOR CC AND NCC CLUSTERS 

The run of k(r) in the handy form Eq. (17) may be inserted 
into the hydrostatic equilibrium Eqs. (8) to yield the Super- 
model for the ICP disposition. Away from the boundary this 
links T(r) and n(r) inversely; such an inverse link is particu- 
larly clear at the very center, as is seen from the dependencies 
T c (k c ) and n c (k c ) obtained from Eqs. (8) on considering the 

explicit scalings T c oc kl ' 5 and n c oc k c 3 ^ 5 with the further fac- 
tor supplied by the dominant integral that scales as £~ a25 (see 
Appendix A). Thus the overall scaling laws read 

T c (xk^ 35 , n c ock- y ; (18) 



FIG. 6. — Same as in Fig. 5 for the CC cluster A2204. The solid line rep- 
resents the Supermodel for a given DM a = 1.27-profile with concentration 
c = 4, and ICP entropy with slope a = 1.1 and central value k c = 10 . The 
dotted line refers to the approximation presented in § 6. 1 . 

these yield T c oc «^ 35 , from which we see clearly how 
low/high T c correspond to high/low n c . In addition, it it is 
easily perceived, and is seen from Fig. 3, that the central runs 
of n(r) are angled or core-like, depending on k c being low or 
high. 

There is much more. Fig. 3 show that the Supermodel 
provides simultaneous, accurate descriptions to both the full 
profiles of T(r) as directly given by Eq. (8), and of the 
surface brightness in X rays provided on integrating along 
the l.o.s. the volume emissivity for optically thin thermal 
bremsstrahlung; this reads 2.4 x 10~ 27 Cx erg s -1 cm -3 with 
jCx °t n 2 T l l 2 (emission lines add for k B T < 2 keV, see Sarazin 
1988). After Eq. (8) one has 

£ x (r)ock-^ Q (r)[l + 2/5b R I \r> r 3 / 5 (r') v1(r')/r'] 7/2 , 

(19) 

with the central scaling given by £x(0) oc £ C L8 . These profiles 
depend: strongly on the value of k c that primarily governs the 
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FIG. 7. — Same as in Fig. 5 for the UNC cluster A1413. The solid line rep- 
resents the Supermodel for a given DM a = 1 .27-profile with concentration 
c = 5.5, and ICP entropy with slope a = 1.1 and central value k c = 2.5 X 10 . 

central pressure and hence the central density run (see Fig. 3); 
weakly on a (and bit) that governs the middle run; mildly on 
the DM concentration c (see § 2) that governs the outer de- 
cline toward the boundary values. 

We illustrate in Figs. 3 and 4 how straightforwardly the Su- 
permodel describes various observables for both classes of 
Cool Core (CC) and Non Cool Core clusters (NCC) as identi- 
fied by Leccardi & Molendi (2008), and also for their interme- 
diate class of UNC clusters. In Figs. 5, 6, and 7 we focus on 
the specific cases of the clusters A2218, A2204, and A1413 
for which both high-resolution XMM-Newton and prelimi- 
nary Suzaku data are available. 

From these results is clearly seen how CC clusters are 
marked out by the presence of a peak of T(r) at r ss 0. 1 -Q.2R 
(or equivalently by T c < 7g). The condition for the peak to oc- 
cur is highlighted on recalling from Eq. (8) that 

T(r)oQk(r)n 2li (r) (20) 

applies, as the ICP counterpart of Eq. (6) for the DM. Given 
that n(r) rises monotonically inward, T(r) will peak and then 



decline toward the center when £ 3 / 5 (r) decreases strongly to- 
ward a low value of k c , as is the case with the CC clusters. On 
the other hand, T(r) will rise to a central, roughly isothermal 
plateau for sufficiently high values of k c . 

From the condition for a maximum to occur in the func- 
tional form of T(r) as given by Eq. (8), on using v 2 c (r) from 
the a-profiles with a = 1.27-1.3 the threshold value for the 
peak reads 

£ e « 2.5 x 10" 2 ; (21) 

with boundary values k.R ~ 1500-2000 keV cm 2 , these cor- 
respond to k c w 40-50 keV cm 2 . In closer detail, the peak 
looms out (the case of UNC clusters) for fc r « 2.5 x 10~ 2 , and 
stands out (the case of CC clusters) for k c < 10~ 2 correspond- 
ing to k c « 15 - 20 keV cm 2 . 

We stress that a finite T c constitutes a natural feature 
of the equilibrium for CC clusters rather than some peculiar- 
ity of cooling flows (see discussion by Peterson & Fabian 
2006). This holds with realistically small values of k c \ but 
even if k c were formally null, T c would decline toward the 
very center with a somewhat flatter powerlaw r 3fl / 5 compared 
to <t 2 oc r 3 "/ 5 . This behavior is consistent with the nonra- 
diative runs of the hydrodynamical simulations by Borgani 
(2007, see his Fig. 1). 

A feature typical of the Supermodel is the peak of T(r) 
closely following the maximum of the DM velocity dispersion 
cr 2 (r), not that of the circular velocity v 2 (r), see Fig. 3 (mid- 
dle right panel). As the former moves considerably down- 
ward with masses ranging from 10 15 to 10 13 M Q (see Fig. 1 
and LC09), we predict the peak of T(r) should also move to 
progressively lower radii in going from rich to poor clusters 
and groups; preliminary data by Nagai et al. (2007) support 
the prediction, but robust data require secure Chandra calibra- 
tions. 

Finally, the Comptonization parameter for the SZ effect 
(Sunyaev & Zel'dovich 1972) obtains from integrating along 
a l.o. s. the volume quantity y oc p in terms of the (thermal- 
ized) electron pressure p = nkgT; after Eq. (8) this means 

y(7) cx [1 +2/5b R J_ dr' r 3 / 5 (r')v 2 (r')/r'] 5 / 2 , (22) 

to the result plotted in Fig. 4. In the central region this scales 
as^(0)cx^ a55 . 

5. STABILITY OF THE CENTRAL CONDITIONS 

Having shown how the Supermodel focuses the central con- 
ditions, here we pursue the discussion ending § 3, and argue 
that they are robust against energy losses or additions. 

5.1. From Cool Cores to cooling cores, and back 

Plainly, the CC state produced by the Supermodel differs 
not only from a cooling flow but also from a freely cool- 
ing core; in fact, cooling is not included in Eqs. (7) as they 
stand. However, the Supermodel focuses the conditions for 
enhanced radiation and fast cooling to set in, namely, low 
though finite T c oc k® 35 linked to high n c so as to imply 

t c « 0.3 ( j] Gyr . (23) 

Vl5keVcm 2 / 

Fig. 3 (bottom left panel) illustrates the run t c (r) of the cooling 
time throughout a CC cluster. 
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So far as it goes for the Supermodel proper. The sequel 
of the story is long accepted in general terms, to imply that 
such an enhanced radiation will lead to entropy loss, that in 
turn will further lower T c and increase n c , so shortening t c 
and opening the way for a classic cooling catastrophe to set 
in (White & Rees 1978; Blanchard et al. 1992). A possible 
happy end to the cooling story in clusters has been widely 
proposed and discussed (see Binney & Tabor 1995; Cavaliere 
et al. 2002; Voit & Donahue 2005; Ciotti & Ostriker 2001, 
Tucker et al. 2007), to the effect that, before the conditions 
run away into a full catastrophe, the ICP condensing around 
central massive galaxies and mixing with their ISM is very 
likely to kindle up AGN activities by renewing mass accre- 
tion onto their powerhouses, the supermassive BHs lurking 
and starving at most galactic centers. Thus recurrent loops 
are conceivably started out by cooling that rekindles AGNs, 
that in turn feed back energy into the surrounding medium 
to eventually quench the BH accretion flow and ultimately 
yield quasi-steady, widespread conditions with k c reset over 
timescales of some 10 8 yr to values around 15 keV cm 2 . 

5.2. Saturation in Non Cool Cores 

On the other hand, NCC conditions prevail above the divide 
at k c ss 40-50 keV cm 2 . We understand the fewer NCC rel- 
ative to the CC clusters primarily on the basis of strong AGN 
inputs, in the tail of the AGN luminosity function N(L) (e.g., 
Richards et al. 2006); this implies for the output statistics 
AEN(AE) oc LN(L) a sharp decline above AE w 5 x 10 61 
erg. We also note that hotter cluster centers tend to im- 
pair or prevent the supersonic condition necessary for driving 
strong pressurized blastwaves (see Appendix B), i.e., M 2 « 
l-AE/E > 3 corresponding to AE /2k B T c m > 1; meanwhile 
the entropy deposited reads k c w AE/{\ - AE/2E) 2 / 3 . This 
will saturate the effects of multiple inputs, if any, occurring 
within a cooling time. 

Similar arguments also apply to the possibly larger if gen- 
erally rarer energy outputs associated to substantial mergers; 
these may be up to AE < 10 64 ergs (see § 3.1.3), but by the 
same token the energies effectively transferred to the ICP in 
NCC conditions are especially prone to saturate after the first 
such event. 

5.3. Concluding on central conditions 

Taking up from § 3.1.3, we recall that the cooling-fueling- 
feedback machinery calls for a correlation between CC clus- 
ters and current central AGNs (the so called dichotomy, see 
Voit 2005) that finds support from a considerable body of ob- 
servations (e.g., Mittal et al. 2009). This picture is currently 
under scrutiny in study cases such as provided by the poor 
cluster AWM4; there the considerable value k c « 60 keV cm 2 
calls for a fueling process unrelated to cooling, with energet- 
ics higher than implied by the current radio activity. In addi- 
tion, no signatures or fossil imprints have yet been found of 
large energy inputs caused over the past 10 _1 Gyr by either 
major AGN outbursts or mergers (Gastaldello et al. 2008, Gi- 
acintucci et al. 2008). This may constitute as of today one 
instance of an exceptional preheating level, standing out as a 
main component to k c . 

To complete the picture, AGNs as widespread agencies for 
rising the central entropies lead to understand the steep de- 
cline of the local L x -T correlation from clusters to groups, 
or the equivalent saturation in groups of the k— T correlation 



(see Ponman et al. 2003; Pratt et al. 2009), in terms of com- 
parable single outbursts in differently sized galaxy systems 
(Cavaliere et al. 2002, LCM05). Finally, AGN kinetic plus 
radiative outputs also lead to expect for clusters a non mono- 
tonic rise and fall of the relation L x - z, consistent with the 
current data (Cavaliere & Lapi 2008). 

In sum, the two main modes for energy injections, namely, 
central AGNs and major mergers, provide different levels of 
entropy input; whence we expect a bimodal distribution for 
the observed number of clusters as a function of the central 
entropy level k c . In fact, the two peaks should be remolded 
to an actual distance considerably smaller than the factor 10 2 
separating the two maximal input levels; this is because the 
statistics will be eroded at low k c by fast cooling, while lim- 
ited at high k c by the small number of strong input events. 
The expected outcome will be not unlike the findings recently 
presented by Voit (2008) and Cavagnolo et al. (2009). 

Back to our main course, we conclude that both the CC and 
the NCC central conditions envisaged by the Supermodel are 
made robust by processes additional, but naturally geared to 
it. 

6. LIMITING MODELS 

While the Supermodel can yield accurate representations of 
the ICP state with a moderate amount of formalism, it is nev- 
ertheless worthwhile to have at hands simple limiting models 
amenable to prompt analytical computations. 

6.1. Mirror dispersions 

We take up the point made in § 4 as to T(r) following a 2 (r) 
for CC clusters, and show in Figs. 3 (middle right) and 6 (bot- 
tom panel) the good performance around the T(r) peak and 
shortward (but for the very central range where T(r) deviates 
upward to its value T c oc k Q c 35 ) provided by the simple model 
with the ICP mirroring the DM dispersion 

T = a 2 1(3 . (24) 

The normalizations are included in the constant parameter 
ft = pm p a 2 /kgT that in Fig. 6 is fixed at 0.75, just the nat- 
ural value it takes when evaluated at R. 

This is similar to the approximation discussed by Cavaliere 
& Fusco-Femiano (1981), and similarly yields the density in 
the explicit form 

n(r) = p t3 (r)a 2W - 1) (r) . (25) 

Here the DM density p{r) is provided by the a-profiles of § 2, 
to imply 

«(r)ocr Q(/3 - 1) /3 5/3/3 " 2/3 ('-) , (26) 

which goes into the simple form n(r) oc p(r) toward the cen- 
ter in the range where p(r) — > r 3 "/ 5 holds but still Eq. (25) 
applies. 

6.2. Polytropic (i-models 

For NCC clusters, instead, we make contact with the classic 
/3-models discussed by Cavaliere & Fusco-Femiano (1978). 

In the central range the contact obtains on noting these clus- 
ters to be marked by high values of k c that cause a nearly flat 
central run of the entropy after Eq. (17). It is now matter of 
straightforward algebra to see that on taking k(r) ss k c to a ze- 
roth approximation, this may be extracted from the integral in 
Eqs. (8), to yield directly 

r/r e = (n/« c ) 2 / 3 = l+2/3 c A0„/5 (27) 
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corresponding to a polytropic approximation with macro- 
scopic index T = 5/3 + dlog A:/dlog n = 5/3. To a next ap- 
proximation T may be obtained by carrying further the ex- 
pansion of the integral, to obtain 



2 1- 
5 ~i 



A0cr Jc dr> 



to be used in the expression 



T/T c = (n/n c y ~ l = 1 + (r- V&Afo/T 



(28) 



(29) 



In the outer regions a similar expansion may be pursued 
both for CC and NCC clusters, to yield an effective index 



5 A(j) r .R J T 



dr' 



(30) 



In fact, Fig. 3 last panel shows that values roughly constant in 
the range T= 1.1-1.2 apply to the outer regions of all cluster 
categories. 5 

Another approximation of the polytropic type is seen to ap- 
ply for r > r m , i.e., to the right of the peak of a 2 (see Fig. 1), 
and may be formally based as follows. Consider that wherever 
density and temperature follow (piecewise) powerlaw runs, 
the elimination of r provides a link of the polytropic form 
nT oc « r ; a similar consideration applies to the DM, leading 
to define a corresponding index A. Thus we may write the first 
and the third sides of Eq. (2) in these terms, and equate them 
directly; then simple algebra provides the explicit relation 



L„ A£^i 
r m r A-i 



(31) 



This, complementarily to Eq. (24), shows that the temperature 
run T(r) tends to follow the DM dispersion cr 2 (r) except for 
the vicinity of the virial boundary where the shock condition 
sustains it at the value T«, and for the very center where a 
finite if small k c matters. 

Summarizing the trend highlighted by the limiting mod- 
els, the passive mirror behavior of the ICP with T oc er 2 pre- 
vails unless is offset by energy inputs, as in fact occurs at the 
boundary for all clusters, and in the central region for the NCC 
clusters. 

7. DISCUSSION AND CONCLUSIONS 

This paper introduces a novel look to the Astrophysics of 
galaxy clusters, in terms of both the a-profiles for the initially 
cold DM and of the Supermodel for the hot ICP. 

As for the DM halos, we have taken up from LC09 the 
physical a-profiles. These are based on the Jeans equilibrium 
between self-gravity and pressure modulated by the DM 'en- 
tropy' run K(r) oc r a . The latter is found from many recent 
Af-body simulations (recalled in § 1) to apply with a closely 
constant within the halo body; we have semianalytically com- 
puted the halo two-stage development and obtained the nar- 
row range a « 1.27-1.3 from poor to rich clusters (see § 2 
for details). The ensuing a-profiles, depending on the two 
key parameters a and c, provide density runs p(r) that sat- 
isfy physical central and outer boundary conditions at vari- 
ance with the empirical NFW formula, and also yield better 

5 For a quick evaluation of the index, consider that V = 5/3 + 
rflog k/d logn = 5/3-a R /g R = b R /g R . 



fits to detailed data from gravitational lensing in and around 
massive clusters (see Lapi & Cavaliere 2009b). 

The ICP, on the other hand, settles to equilibrium within 
the gravitational wells associated with the a-profiles, under 
control from the thermodynamic entropy produced by bound- 
ary and central shocks driven by AGNs or major mergers, 
plus a possible preheating basal level (see § 3). These phys- 
ical effects may be compounded in the two-parameter form 
k(r) = k c +(l—k c )r a with a ranging from 0.8 to 1.1. The result- 
ing equilibrium may be concisely rendered as a trend for the 
ICP to follow the DM in the passive behavior T(r) oc a 2 (r), in 
the radial range free from the energy inputs that steadily pro- 
duce the outer boundary slope and intermittently refresh the 
central level k c . 

In detail, our Supermodel of § 4 provides accurate and ex- 
tended representations for the runs of ICP temperature and 
density and of the related ICP observables (see Fig. 5, 6, and 
7 for examples). These validate the assumption of hydrostatic 
equilibrium, and closely constrain not only the value of the 
concentration c for the DM a-profile but also the two ICP 
parameters a and k c . In fact, such representations hold for 
either Non Cool Core and Cool Core clusters, marked out, 
respectively, by a monotonic outer decline of T(r) from a cen- 
tral plateau at T c > T«, or by a middle peak. In the Super- 
model these morphologies are produced by the central entropy 
k c being higher or lower than a threshold value k c « 20 - 50 
keV cm 2 ; correspondingly, the X-ray brightness features a flat 
corelike or a steep central run, based upon the same mildly 
cusped DM a-profile. In a forthcoming paper, we will present 
a detailed analysis of an extended sample of NCC and UNC 
clusters, with the aims of disentangling the origin of the cen- 
tral energy inputs, and of evaluating the variance in the cluster 
ages through their outer DM concentrations (see § 2). 

We stress that the Supermodel links T(r) and n{r) inversely 
from the bulk toward the central region (see § 4); in particular, 
it yields for the very central values the scaling T c oc fcj! 35 and 
n c oc k' 1 , implying t c oc k l c 2 for the cooling time. The stability 
of such values is argued in § 5. High T c combines with flat 
n c to produce in NCC clusters central conditions conducive to 
saturation of k c toward values around 10 2 keV cm 2 . In CC 
clusters, instead, low though finite T c combines with high n c 
into a condition paving the way to fast cooling; this condition 
is conducive to triggering intermittent, recurrent loops going 
through the stages: cooling, massive BH fueling, AGN energy 
feedback, that halts further fueling and activity; these loops 
make possible in the long term a quasi-steady state. In the 
center of NCC clusters, instead, hot conditions suppress AGN 
reactivations owing to lack of dense cool ICP crowded around 
the central cluster galaxies; in addition, they tend to saturate 
the effective energy coupling from the most powerful AGNs 
or mergers by preventing or impairing supersonic conditions 
conducive to strong shocks. 

Concerning central conditions, we emphasize two points. 
First, in CC clusters the Supermodel predicts a finite (non- 
zero) central T c with no need for any twist in cooling flow 
theories; this rather constitutes a natural condition set by their 
equilibria at low k c and stabilized against cooling by recurrent 
AGN activity. Second, in NCC clusters we expect saturation 
to enforce stability of the higher k c levels set by the inputs 
from powerful AGNs and from the stronger if rarer mergers. 

Moving into the middle radial range where energy sources 
may be neglected, the 'mirror' behavior of Eqs. (24) and (31) 
prevails with T(r) passively following cr 2 (r); a novel feature 
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FIG. 8. — Reconstruction of the entropy profiles from joint X-ray and SZ 
observations (Cavaliere et al. 2005, see their Eq. [8]). The dashed line shows 
the DM entropy profile with the slope a = 1.3. Solid lines are the input ICP 
entropy profiles described by Eq. (17), for our standard values k c = 5 X IQ~- 
(red), k c = 2.5 X 10~ 2 (cyan), and k c = (blue); the error bars illustrate the 
reliability of the reconstruction from joint mock observations in a rich cluster 
at 1 Gpc of the X-ray brightness with resolution of 1 " and sensitivities AIx = 
2 X lCT 4 Ixc, and of the SZ effect with resolution of 10" and sensitivities 
A/ v = 10~ 2 The blue dotted line represents an alternative ICP entropy 
profile with a = 0.8 (and k c = 0) as we expect for ICP preheated at 1/2 keV 
per particle (see § 3.1.1); similar shapes have been observed by Sun et al. 
(2009) in many poor clusters, and by Lemze et al. (2008) in A1689. 



emerging from the Supermodel is the very close location of 
the two respective peaks. This is because such two homolo- 
gous quantities arise from a parallel response to the require- 
ment of withstanding the common gravity for equilibrium 
(see § 6, also Fig. 3). This gravity-induced behavior is at 
the root of the remarkable effectiveness of the simple model 
T(r) oc cr 2 (r), which for the CC clusters holds well down to- 
ward the center (see Fig. 6). As a consequence, the Super- 
model predicts the peak of T(r) to move toward progressively 
smaller radii in going from rich clusters to groups. 

On approaching the boundary, the run of T(r) again de- 
viates upward from this passive trend (see § 3.1 and Ap- 
pendix B), with a boundary value Tr sustained by the en- 
ergy input associated to infall. Here the passive ICP behav- 
ior is again expected to be broken by the energy transfer 
due to electromagnetic interactions and localized to a range 
Ar ~ X pp <C r (whilst bubbles or shocks starting from the 
center smear their energy out to some 10 2 kpc, implying an 
effective Ar/r ~ 1). However, pinning down the outer devia- 
tions requires high sensitivity, currently achievable only with 
full use of the Suzaku capabilities. 

Next we highlight an unexpected connection specifically 
emerging from the Supermodel. X-ray observations of clus- 
ters yield information concerning the concentration c through 
the values of the outer entropy slope a, and more directly from 
detailed fits to the surface brightness data; in fact, we expect a 
to be lower and the density profile to be steeper for early clus- 
ters with higher concentrations c (see § 3, in particular below 
Eqs. [13]). This specific prediction may be tested through 
extended simulations covering high-c halos and nonadiabatic 
processes as discussed, e.g., by Borgani (2007). In paral- 
lel, high concentrations are keenly pinpointed by gravitational 
lensing observations (Lapi & Cavaliere 2009b and references 
therein). This opens the way to the use of existing X-ray data 
as convenient pointers to targets for time-expensive gravita- 



tional lensing observations. 

To conclude, we turn to contrast the ICP and the DM behav- 
iors. Note that the two density runs p(r) oc [a 2 (r) /K(r)] 3 ^ 2 and 
n(r) oc [kBT(r)/k(r)] i ^ 2 will differ even where T oca 2 applies, 
to the extent that the DM and ICP entropy runs differ. This 
brings us to directly compare these two governing entropies. 

In a nutshell, their common features stem from smooth, 
slow gravitational mass infall onto the outskirts, while their 
detailed runs both in the outer and in the central range re- 
flect their different sensitivity to other energy inputs. Quan- 
titatively, both the underlying key parameters bu and K cr i t 
are amenable to conversion of infall kinetic energy. They 
take on very close values w 2.65-2.55 (see § 3) and 
ft cr it ~ 2.6-2.5 (see § 2) at their corresponding fiducial points 
r « R or r = r p ; for increasing a or concentration c, they de- 
crease together since both depend on 1 / Atfi in terms of the 
relevant potential drops from the turning point to R or r p (see 
Eq. [10] for the ICP and Eq. [13] in LC09 for the DM). 

On the other hand, the differing features stem from local 
vs. non-local character of the energy conversion. The colli- 
sionless DM particles fall from the cluster surroundings well 
into the body, where their kinetic energy is non-locally and 
progressively randomized, and spreads out entropy by orbit 
superposition and stratification with widely distributed apoc- 
enters (see LC09 and references therein). Correspondingly, in 
the DM halos K(r) = a 2 /p 2 / 3 starts out in the outskirts with 
uniform values a ~ 1.1 in all clusters, to steepen toward the 
body to universal values 1 .27- 1 .3 in a gently convex shape 6 . 

Meanwhile, in the ICP k{r) = k B T/n 2 / 3 stalls out at the 
boundary r s=s R from somewhat lower average values a w 
1.1 — 1.2, but with a considerable variance when large pre- 
heating and high DM concentrations are included, see § 3.1.1 
and evidence referred to therein. The slope, if anything, flat- 
tens out toward an effective value a m in the presence of 
central energy inputs. 

Beyond details, this concave (vs. convex) shape of a toward 
the center, together with the sensitivity of its bondary values 
to outer potential and preheating constitute features specific to 
the ICP. Thus we conclude that basically similar gravitational 
processes in DM and ICP (randomization of bulk kinetic en- 
ergy, but on different scales) with the addition of the ICP col- 
lisional sensitivity to other energy inputs, concur to produce 
dissimilar shapes for K(r) and k(r). 

Model independently, we propose two observational tests 
addressed at directly probing in clusters the two underlying 
entropies. The run K(r) of the DM entropy can be derived 
from probing the a-profile by gravitational lensing observa- 
tions as recalled in § 2. The run k(r) of the ICP entropy can 
be reconstructed as proposed by Cavaliere et al. (2005) and 
illustrated in Fig. 8 starting from the relation 

k(r) = y^(r)£- x 10/9 (r), (32) 

that joins deconvolved observations of X-ray brightness and 
SZ effect irrespective of any modeling or assumption on hy- 
drostatic equilibrium and of redshift information (Cavaliere 
et al. 1977). Such studies may be particularly useful in the 
ongoing search for early clusters (see Andreon et al. 2009). 

6 Note that our slope of K = a 2 / p 2|/3 defined in terms of the 1-D velocity 
dispersion is consistent with that in terms of the 3-D dispersion, e.g. Fal- 
tenbacher et al. (2007); in detail, Dehnen & McLaughlin (2005) and Ascasi- 
bar & Gottlober (2008) find the latter to be somewhat flatter than the former 
in the outer body where radial anisotropics tend to prevail. 
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TABLE 1 
Fit parameters of Eq. (Al) 



kc 




c = 3.5 






c = 4.5 






c = 5.5 




Ao 


Ai 


A 2 


Ao 


A, 


A 2 


Ao 


Ai 


A 2 





17.77 


2.86 


0.49 


22.82 


2.72 


0.47 


27.59 


2.64 


0.47 




17.47 


2.86 


0.49 


22.33 


2.71 


0.48 


26.88 


2.63 


0.48 


l(r 3 


15.93 


2.83 


0.54 


19.99 


2.66 


0.54 


23.72 


2.56 


0.53 


5 x ur 3 


13.42 


2.73 


0.62 


16.53 


2.51 


0.61 


19.34 


2.38 


0.61 


ur 2 


12.04 


2.63 


0.66 


14.70 


2.39 


0.66 


17.09 


2.25 


0.65 


2.5 X l(r 2 


10.06 


2.44 


0.72 


12.14 


2.18 


0.71 


13.99 


2.02 


0.71 


5 x ur 2 


8.53 


2.26 


0.77 


10.19 


1.99 


0.75 


11.66 


1.83 


0.74 


7.5 X 1CT 2 


7.64 


2.14 


0.79 


9.08 


1.87 


0.77 


10.34 


1.71 


0.76 


10-' 


7.03 


2.05 


0.81 


8.31 


1.78 


0.79 


9.44 


1.62 


0.77 



NOTE. — In computing the integral Eq. (8) we have used the DM a-profiles with a = 1.27. The above fitting coefficients are given for three values of the 
concentration parameter; other values may be derived by standard interpolation techniques. 
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APPENDIX 
THE INTEGRAL IN EQ. (8) 

Here we provide an effective approximation to the integral I(r) appearing on the right-hand side of Eq 
following analytical expression 

I(r) = \ %■ v 2 c {?) r 3 / 5 (r') ~ A exp(-A l r A ^)- 

J r ^ 

the fitting parameters Aq,A\, andA2 depend weakly on k c , as specified in Table Al. The dotted lines in Fig. 2 (top panel) illustrate 
the effectiveness of such an approximation for our three standard values of k c . Note that on approaching the center the integral is 

-1/4 

numerically found to scale with central entropy k c like I c cck c , see Fig. 2 (bottom panel). This result is used in § 4 of the main 
text. 

We plan to provide elsewhere a set of fitting formulae for p(r) and v 2 (r) leading to an analytic expression for I(r), as tools 
enabling direct data analysis and extensive precision fits. 



. (8), in terms of the 
(Al) 



ENTROPY PRODUCTION IN SHOCKS 



While cooling may condense out the colder fractions of the ICP and indirectly raise the average entropy of the rest, it is 
generally agreed (see Cavaliere et al. 2002; Voit 2005) that substantial entropy production requires Shockwaves. Here we recall 
from LCM05 the temperature, density and entropy jumps produced across a shock transitional layer. 

Conservation of mass, energy and total stress across the latter lead to the classic Rankine-Hugoniot temperature jump 



— = — M +-■ 
Ti 16 8 



3 1 
16 M 2 



(Bl) 



The subscripts 1 and 2 denote the pre- and post-shock quantities, and M = (3 \im p v\l 'SkaTi) 1 ' 2 the shock Mach number; the 
quantities with tildes refer to the shock reference frame, which is convenient in the case of internal shocks driven, e.g., by AGNs. 

On the other hand, in the case of accretion shocks it is more convenient to work in terms of the infall velocity v\ and of the 
related Mach number Al = (3 i±m p v\/5kBT\) x > 2 measured in the center of mass frame. Assuming the downstream kinetic energy 
to be small, one finds the temperature jump in the form 



72 

7i 



4 , 
1 + -M 2 



1 



1 + - 



1 

4 M 2 



In either reference frame the density jump in terms of the pre- and post-shock temperatures reads 




(B2) 



(B3) 



The entropy jump K 2 /K\ = (T 2 /Ti)/(n 2 /ni) 2 ^ may be easily composed from the relations above. 
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Handy expressions apply to strong shocks, when the above expressions reduce to 

T 2 3 nm p v\ 7 1 (xm p v 2 3 Hi ^ 4 ( ■ (B4) 
Ti 16 k B Ti 8 3 k B Ti 2 ' n x \ 16 T 2 ) ' V 

these approximations actually apply to high/intermediate M 2 > 3; the corresponding entropy jumps read 

f<2 3_ _\_ fJMp v\ 3 _1_ 1 _1_ /imp v, 17 _1_ 

^" 16 4273 ~^7T + 2 4^ ~ 3 4273 ~^7r + Y 4273 ' ( ' 

In last relation, the second term on the r.h.s. expresses the contributions of the advected external entropy 0.&4K\, and is relevant 
when k B T\ > Q,\6m p v\, that is, when relatively strong preheating affects the gas infalling into a poor cluster or a group. 

When a central energy pulse is discharged into the equilibrium ICP with a density gradient, e.g. n(r) oc r" 2 , a blast wave is sent 
out, that is, a non-linear perturbed flow that terminates into a leading shock and comprises bulk kinetic energy up to matching 
the thermal one (see Sedov 1959). When the pulse is short compared with the transit time and the effects of gravity and upstream 
pressure are neglected, the Mach number decreases radially as Ai(r) oc r~ l l 2 . When the pulse is sustained during the transit 
time, also A4{r) will be. Even with gravity and upstream pressure considered, this is found to hold in the simple case of a pulse 
constant over the transit time through the region where n(r) oc r~ 2 applies, and to hold also for other combinations of equilibrium 
gradients and pulse shapes (LCM05). For longer times, A4{r) declines and the blast dissipates its kinetic energy into the ICR 

In such cases the Mach number depends on the energy input AE relative to the ICP total energy E in the affected volume, 
simply as M 2 w 1 + AE /E; the condition for a strong shock M 2 > 3 clearly translates into AE /2mk B T > 1. It turns out that the 
ICP may be (partly) evacuated from the central region to a residual average density n(l — An/ri) oc 1—AE/2E, with an associated 
entropy k c oc AE/(\ - AE/2E) 2 ' 3 . 
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